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I. INTRODUCTION 



Thanks to recent development of technology, it has become possible to 
synthesize many quasi-low-dimensional systems experimentally. This has 
led to growing theoretical interest in low-temperature properties of low- 
dimensional systems. A remarkable example is the Haldane conjecture 
Hill for the spin-S 1 antiferromagnetic Heisenberg chain. 

Recently, two of the present authors considered |3| low-temperature 
low-magnetic-field properties of the spin-S 1 ferromagnetic Heisenberg 
chain. The Hamiltonian of the model is defined by 

J L h L 

ft = ~-^J2 s i ■ s i+i ~ t;H^ 2 ' 
D i=i D i=i 

where the coupling constant J is positive and h denotes a magnetic field. 
They calculated |J the free energy and the magnetization of the S — 1/2 
system numerically, employing the thermal Bethe-ansatz method |4| for 
the thermodynamic limit L —>■ oo. They found || that the magnetization 
for the 5 = 1/2 system has the same scaling function as that of the 
classical (S = oo) ferromagnetic Heisenberg chain: 

m(T, h) ~ m (x) , (1.2) 

2 44 / \ 

= -x - — x 3 + O (x 5 ) . (1.3) 

Here x is the scaling variable: 

Jh , . 

They thereby conjectured that this scaling function should be universal, 
or common to all values of S. This conjecture implies that the spin- 
wave excitation from the ground state of the ferromagnet show universal 
behavior. 

In order to confirm this conjecture for 1 < S < oo, we have to obtain 
information from numerical calculations of finite systems, because the 
Bethe ansatz method is not applicable to higher-spin Heisenberg chains. 
In the previous paper, we hence generalized the above scaling form 
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( |1 . 2| ) to finite-size scaling and conjectured that the scaling function of 
the magnetization, 



m(T, h, L) ~ m (x, y) 



;i.5) 



is universal for the arbitrary-S 1 ferromagnetic Heisenberg chain. Here the 
scaling parameters x and y are defined by (|1.4|) and 



y 



J 

TL' 



1.6) 



It is difficult to obtain the magnetization m(T, h, L) for the classical 
system. Hence, as the first step, we treated || the finite-size scaling 
form of the linear susceptibility. The scaling form ( |1.5|) is followed by the 
finite-size scaling function of the linear susceptibility Xi in the form 



dm 



h=0 



J _ / J 



(1.7) 



where 



xi (y) 



drh(x, y) 



dx 



In the previous paper, we analytically obtained the function Xi f° r 
the classical, or the S = oo Heisenberg chain with both periodic and 
open boundary conditions. Moreover, we showed that the function xi 
for S = oo fits numerical data for 5 = 1/2 and S = 1 quite well. Thus 
we partly confirmed the conjecture on the universality of the scaling 
function (|1.5| ). 

In the present paper, we further treat the scaling function of the 
nonlinear susceptibility, and show its universality. Equation (\1.3j) ensures 
the universality not only of the linear susceptibility xi but also of the 
third-order nonlinear susceptibility Xz'- 



1 ^ 



m 



3! dh 3 



h=0 



J 3 . /J_ 

Te X3 y TL 



1.9) 



where 
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1 d 3 m(x, y) 

Xs \y) ~~ 



(1.10) 

x=0 



3! <9x 3 

The result (|1.3|) in the thermodynamic limit requires the relation 
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X 3 (l/ = 0) = - — . (1.11) 

In the present paper we obtain an analytical expression of %3 f° r the 
classical chain. (This was partially reported in the previous paper. ||) 
We then show that the function X3 obtained for 5 = oo is consistent with 
numerical data for quantum systems, namely, 5 = 1/2 and 5 = 1. Thus 
we confirm the universality of the scaling function (|1.5| ) up to the third 
order of x. 

In the course of the study, we also obtain the exact expressions of the 
four-point correlations and the nonlinear susceptibility of the classical 
system. To our knowledge, these expressions for the periodic classical 
system have not appeared in the literature. 

Before going into details of calculations, we comment here on the 
definition of the scaling limit. We see in the following that there are 
corrections of the form yy/T, yT, etc. to the finite-size scaling functions 
(1.7) and (|1.9|) . These are less singular than the leading term of y for 



large systems at low temperatures. We have to exclude these correction 
terms in order to obtain the scaling function of y. We hence extrapolate 
the limit T — > while fixing the scaling variable y = J/ (LT) . We 
refer to this limit as the scaling limit. It should not be confused with the 
thermodynamic limit y — > 0, or the limit L — > oo with T fixed. Even after 
taking the scaling limit, we still consider a finite value of y = J /(TV). 

We organize the paper as follows. In §2 we analytically calculate the 
finite-size scaling function of X3 of the classical Heisenberg model for the 
periodic and open boundary conditions. We show numerical confirmation 
of the universality of the scaling function for 5=1/2 and 5 = 1 in §3. In 
§4 we discuss origins of the universality of the third-order nonlinear sus- 
ceptibility. In Appendix A we present in detail the analytic calculation of 
the four-point correlation function of the classical Heisenberg chain with 
periodic boundary condition. In Appendix B we show the derivation of 
the third-order nonlinear susceptibility of the classical Heisenberg chain. 
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II. FINITE-SIZE SCALING FUNCTIONS OF THE 
THIRD-ORDER NONLINEAR SUSCEPTIBILITY 



In the previous paper || we presented the finite-size scaling func- 
tion of the linear susceptibility, xi- We used the classical Heisenberg 
chain to obtain the analytic form of Xi- m this section we derive the 
finite-size scaling function of the third-order nonlinear susceptibility, xzi 
analytically for the classical Heisenberg chain with the periodic and the 
open boundary conditions. We insist that the ferromagnetic Heisenberg 
chain has this finite-size scaling function X3 not only in the classical case 
(S = oo) but also in the quantum case. 

The classical Heisenberg chain is given by 

J L h L 

D 1=1 D i=\ 

L L 

— > — J rij • THfi — h n\ as S — >• oo, (2.1) 

i=i i=i 

where { } are vectors of length unity. The nonlinear susceptibility is 
given by summation of four-point correlations as follows: 

Xa(T,L) = 31^3^4 E [( S ii S l S i 3 S l 

h,i2,i3M 
cz qz \ / qz qz 



Q2 \ / QZ QZ 



/ QZ QZ \ / QZ QZ 



(2.2) 



where (■ • •) denotes the thermal average. 



A. Periodic boundary condition 

First, we consider the periodic system. To obtain the nonlinear sus- 
ceptibility, we have to calculate the partition function, the two-point 
correlation function and the four-point correlation function. 
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Joyce || obtained the partition function of the periodic chain, Z^, in 
the form 



Z L {K)=Y J {U + l)\ l {K) 1 
1=0 



;2.3) 



where K = J/T and 



7T 



(2.4) 



with I l+ i(K) being modified Bessel functions of the first kind. The two- 
point correlation function schematically shown in Fig. |l| is given by HH 



z z \peri 

n'h+M/L 



3Z 



L 1=0 



u 



M 



U 



L-M 
I 



(2.5) 



where the bracket (• ■ ■)^ n denotes the thermal average for the periodic 
system of length L, and the function ui is defined by 



ui{K) 



(2.6) 



From the calculations given in Appendix A, we also obtain the four-point 
correlation function shown in Fig. |^ as 

j z z z z \peri 

\ n l n l+A n l+A+B n l+A+B+CIL 
i oo 



J L 1=0 



ufvfuf + ufvfu? + ufvfuf + ufvfuf 



u?u D 



} 



+9i [ u f u ? + ■«*! ■«*» 
where D = L- A- B -C, 

2 (/ + !)(/ + 2) 



fx 
9i 



15 11 + 3 
1 (/ + l)(4/ 2 + 8/ + 5) 
15 (2/ + l)(2/ + 3) 



(2.7) 

(2.8) 
(2.9) 



and 
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vi(K) 



Xi+2(K) 



(2.10) 



These expressions ( |2.3| )-( [27T0| ) give the exact expression of the non- 
linear susceptibility Q2.2|) . See Appendix B for details. Now we calculate 
its scaling form. 

In the scaling limit, that is, in the limit T — > with y = J /(TV) 
fixed, the partition function is expressed in terms of y as follows: 0] 



Zl(K) 



oo / \ 



\o(K) L 
-+W( y )=£(2/ + l)exp 



1=0 

oo 



1=0 



1(1+1) 

2y 



(2.11) 



Next, we can obtain the finite-size scaling function ( |1.9| ) of the third-order 
susceptibility with the periodic boundary condition, xlT^ as follows: 

' X r\T,L)T°- 



xT\y) 



lim 

T-»0 
y fixed - 



J 3 



— y^exp 

my) to 



i(i + i) 
2y 



8 1 + 2 
15 21 + 3 



x 



+ 



(I + 2)y 
1 



exp 



exp 



x 



(2/ + 3)(/ + 1) 1 

Z + l 
(2/ + 3)(/ + 2) 2 

2 4/ 2 + 8/ + 5 
15(2/ + 1) (21 + 3) 

^1 + exp 
1 — exp 



.1+1 

y , 

— exp 



1+1 

y , 

21 + 3\ 



l + l 



y 



exp 



[(I + 1)2/ 
2 / 



/ + ! 

2/ , 



(1 + I) 2 v 
1 2 

y%) 29 ?/' 



(2.12) 
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The key in the calculations of the above scaling forms is to note 
the behavior of the functions u\ and v\. These functions behave at low 
temperatures as follows: M 



ui{K) 
vi(K) 



1 - 
1 - 



l + l 

K 
2/ + 3 

K 



(2.13) 
(2.14) 



Hence they converge to unity as T — > 0, or K = J/T — ► oo except for 
the terms u\ and vf. Using the scaling parameter y = K/L, we have the 
following asymptotic forms of uf and vf in the scaling limit: || 



ui(KY 



vi(K) 1 



1 - 
exp | 
1 - 

exp 



I + 11 

z + r 

y 2 

21 + 3 1 



[2.15) 



2/ + 3 
1/ 



L + F 



12.16) 



Further details are given in Appendix B. 

In the thermodynamic limit y — > 0, the right-hand side of eq. (|2.12|) 
is reduced to 



xT\y - o) 



44 



135 



(2.17) 



This value —44/135 is consistent with the value in the thermodynamic 
limit, eq. Ql.llj) . 



B. Open boundary condition 



Next we consider the open system of length L. Because we handle 
only three modes I = 0, 1, 2 of the modified Bessel functions in the open 
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systems, it is easier to calculate the correlation functions in the open 
systems than in the periodic systems. The partition function is identical 
to unity. Fisher gave the two-point function as follows: 



n'n+M 



open 



-M K ) 



M 



(2.18) 



Tomita and Mashiyama H gave the four-point function as follows: 



( n l n l+A n l+A+B n l+A+B+c)° PCn — g U o(^y 



5 



u (Kf. 
(2.19) 



Using these expressions, we already reported H the finite-size scaling 
function of the third-order susceptibility with the open boundary condi- 
tion for the classical Heisenberg chain as follows: 



xT\y) 



44 32 



exp 



135 45 
169 + 43 exp 



y, 



2y_ 
405 



1 — exp 



y, 



2 exp — 



[2.20) 



In the thermodynamic limit y — > 0, we have the value of the finite-size 
scaling function of the third-order susceptibility as 



xT a (y - o) 



44 



135 



(2.21) 



This value is consistent with the eq. ( |1 . 1 1| ) as well as in the periodic case 

(ra- 



in. FINITE-SIZE CALCULATION FOR 5 = 1/2 AND S = 1 

In the previous section, we obtained the finite-size scaling function 
X3 for the classical systems with the periodic and open boundary con- 
ditions. We conjecture that this scaling function is common to the 
ferromagnetic Heisenberg chain with an arbitrary magnitude of the spin. 
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In this section, in order to check this conjecture, we calculate the third- 
order susceptibility ( j2.2|) of finite systems with 5=1/2 and S = 1 by 
the Householder method and compare the numerical data in the scaling 
limit with the finite-size scaling function (|2.12|) and (|2.20|). 



A. Periodic boundary condition 

First, we treat the periodic system. We diagonalized the Heisenberg 
Hamiltonian (|1.1|) numerically by the Householder method. The system 
size is up to L = 14 for S = 1/2 and up to L = 10 for 5 = 1. We 
then calculated the spin correlation functions and summed them up over 
all sites following eq. ( |2.2| ). We thus obtained the third-order nonlinear 
susceptibility numerically. 

We plot the numerical data of S = 1/2 chains with the periodic 
boundary condition in Fig. [| and those of S = 1 chain with the periodic 
boundary condition in Fig. [|. We observe corrections to finite-size scaling 
in these figures. In order to take the scaling limit, or to exclude the 
corrections, we here assume that the leading correction is of the same 
form as in the linear susceptibility: 



j 

Xs{T,L) = — 



(3-1) 



Uv) + o(y/¥/j) 

In the previous paper B we take the scaling limit, fitting the data 



quadratically with respect to JT/J. In the present analysis, however, 
the quadratic fitting resulted in spurious scaling limit. Some of the ex- 
trapolation curves for the quadratic fitting have the maximum around 
Jrjj ~ 0.2 and start decreasing as T/J — > 0. This indicates that 
corrections higher than the second order are large in the nonlinear sus- 
ceptibility. We therefore used corrections to finite-size scaling up to the 
third-order term, i.e. (T/J) 3 / 2 . Each fitting curve in Figs. |3] and || was 
determined by the four points nearest to the ordinate. The crosses on 
the ordinate of Figs. |3| and |] denote the scaling limit Xs er \y)- 

In Fig. [|, we summarize the data in the scaling limit obtained in Figs. 
|3] and |] together with the finite-size scaling function for the classical case 
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( P-12Q . (In calculating eq. ( |2.12| ) we numerically summed the terms up 
to I = 100, which turned out to be sufficient for the convergence of the 
series.) The data for S = 1/2 and S = 1 are quite consistent with the 
finite-size scaling function for S = oo. This suggests that the finite-size 
scaling function ( p. 12 ), which is given analytically for the classical chain, 
does not depend on the magnitude of the spin S. 



B. Open boundary condition 

Next, we treat the open chain. We perform in Figs. |6| and |7| the same 
analysis as in the periodic systems, i.e. Figs. [| and f|. The finite-size 
scaling function ( [2.20|) and the scaling-limit data from Figs. || and |7| are 
summarized in Fig. We conclude that the finite-size scaling function 
( p.20| ) is common to all magnitudes of the spin for the open systems as 
well as for the periodic systems. 



IV. DISCUSSION 

By using the numerical calculations for the quantum case, namely 
S — 1/2 and S — 1, we showed that the finite-size scaling function % 3 
derived from the classical case is universal with respect to the magnitude 
of the spin both for the periodic and the open boundary conditions. 
The scaling function % 3 is exactly given by eq. ( |2.12j ) and eq. ( |2.20| ) 



for the periodic and the open chain, respectively. We emphasize that 
the universality both of the linear and the third-order susceptibilities 
implies the universality (|1.5| ) of the magnetization of the arbitrary-spin 
ferromagnetic Heisenberg chain. 

The universality of the scaling function implies that the spin-wave 
excitation in the present systems show universal behavior. The ground 
state of the arbitrary-S 1 Heisenberg ferromagnet in a magnetic field is 
obviously the state where all the spins are aligned in the direction of the 
field. The thermal fluctuation creates spin-wave excitations. The present 
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universality indicates that the thermally excited spin wave behaves in- 
dependently of S. We can see this universal behavior as the universality 
of the correlation functions. We have already suggested || that the uni- 
versality of the linear susceptibility originates in the universality of the 
two-point correlation function. We naturally expect that the universal- 
ity of the nonlinear susceptibility originates in the universality of the 
two-point and the four-point correlation functions. 

Indeed, we can derive the scaling forms of the correlation functions 
as follows. The finite-size scaling functions for the linear susceptibility 
and the nonlinear susceptibility are rederived from these scaling forms 
below. When the distance between two spins is large, we have the scaling 
function of the two-point correlation function for the periodic chain as 
follows: § 




pen 



£g{e + i)«p 



exp 




(4.1) 



for W S and M, (L — M) ^> 1, where the correlation length 6 of the mode 
/ is defined by 



6 



K 

ITT' 



(4.2) 



We also obtain from eq. ( p.7|) the scaling from of the four-point correlation 
function for the periodic chain as follows: 




s j \ s j v s 



— y~! ex p 



x {fl 




2y 



pen 



exp 



exp 



l +A+B+ C 
S 



B 



exp 



6+i 



exp 




12 



+ exp 
+ exp 
+ exp 



iy 

6, 



exp 



exp 



C 



■— I ex P I - t- I ex P 



+ <7z 



exp 



6, 



exp 



6+1, 
6+1 



exp 



exp — — 



exp 



.4 



IT 

6, 



exp I — — I exp 



— I + <'xp I ) exp f — — 



6 



a, 



6 



(4.3) 



for V S and A, £, C, D > 1 with A + B + C + D = L. 

For the open chain, we have obtained the scaling function of the two- 
point correlation function as [H] 



Sf 

s 



's 



open 



l+M 
S 



-exp(-M/f ) 



(4.4) 



for V S' and M ^> 1. We may also obtain the scaling function of the 
four-point correlation function as 



Si 

s 
1 



's 



1+A 



's 



l+A+B 



open 



■A+B+C 



S 



S 



s 



f A 

~ - exp — — 
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^0, 



B 



exp I — I exp 



-ir +1 



exp rs, 



(4.5) 



for V S and A 5, C> 1. 
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Appendix A:. FOUR-POINT CORRELATION FUNCTION OF 
THE CLASSICAL HEISENBERG CHAIN WITH THE 
PERIODIC BOUNDARY CONDITION 

In this Appendix we show the derivation of the four-point correlation 
function Q . 

We consider the classical Heisenberg chain ( [2.1|) with the periodic 
boundary condition. We choose the direction of the magnetic field 
h(|h| = h) as the z axis of the three-dimensional spin space. We de- 
fine the polar coordinate 0, and of the ith classical spin relative to the 
z axis. We express each spin in the form 

nl = cos 6i 




(A.l) 



where Y\ m is the spherical harmonics. The exponential function of nj-rij+i 
is expanded in terms of the spherical harmonics as follows: || 

e Knin i+1 _ e /fcos0 

oo I 

= 4vrE E \iYUWU* + l), (A.2) 

1=0 m=-l 

where the parameter B is the angle between the spins n^and n i+ i, and 
A; is defined by eq. (|2.4| ). 

In the following algebra, we use two formulas concerning the integra- 
tion of the spherical harmonics. First, the integral of the product of two 
spherical harmonics over the solid angle Qi is given by |J 



14 



2tt rn 



JO 



(A.3) 



where 5 is the Kronecker symbol. Next, the integral of the product of 
three spherical harmonics can be expressed in terms of Wigner's 3n — j 
symbols: MM 



2- 



X 



-1) 
/ 

V 



(2/ 1 + l)(2/ 2 + l)(2/ 3 + l) 



^i h h 




47T 

n h h 



1/2 



—mi m 2 m 3 



(A.4) 



where Wigner's 3n — j symbols can be written as follows: 
^ h h h ^ 



mi m 2 m 3 
= 6. 



J mi+m2+m3,0 



(h + l2-h)Kh-k + h)K-h + k + hy. 



(h + k + h + iy. 

x [(/1 + mi )\(h - mi)!(J a + m 2 )!(/ 2 - m 2 )!(Z 3 + m 3 )!(Z 3 - m 3 )!]" 



X 



z\(h +I2-I3- z)\{h - mi - z)!(Z 2 + m 2 - z)!(Z 3 - Z 2 + mi + z)\(l 3 - h - m 2 + z)\ ' 

(A.5) 



The summation in eq. ( |A.5|) is over all integers z. 

Using eqs. ( |A.l ) - ( A.4j) , we can write down the four-point correlation 
function shown in Fig. ^| as 

F(A, B,C,D:K) = Z L {n{n\ +A n{ +A+B n^ A+B+c ) l , (A.6) 



2 L 

n 



An) f = V 47r 



da 



exp(irn J n i+1 )F 10 (l)Fi (l + A) 
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xF 10 (l + A + B)Y W (1 + A + B + C) 

3 ^ 2 



J dfii J dVt 1+A J AVt 1+A+B J dtt 1+A+B+ c 



An 

I! I! I! I! K^h^u 

{hmi} {hm 2 } {hm 3 } {I4IJ14} 

x^ m4 (l)y 10 (l)^ mi (l) x Y* mi (l + A)Y l0 (l + A)Y hm2 (l + A) 

xY* m2 (1 + A + B)Y 10 (1 + A + B)Y hm3 (1 + A + B) 

xY* m3 (l + A + B + C)y 10 (l + A + B + C)Y km4 (l + A + B + C) 

3 ) 2 e e e e kkkk 



Att 
x(-l) 



{iimi} {hm2} {/3"i3} {/4"i4} 

mi+m2+m3+m4 ggfi + l)(2/ 2 + l)(2/ 3 + 1)(2Z 4 + 1) 



x 



X 



10 



z 4 



^2 1 / 3 





Z 4 1 h 

— m 4 mi 

k 1 / 3 

— ffl2 1713 



hi 




h 1 ^2 

—mi m2 



'3 1 <*4 





t 3 1 u 

-m3 m 4 



(A.7) 



In order to calculate eq. (|A.7|) further, we need the following formula for 
the Wigner symbol: M 



(A.8) 







( 


I 11'^ 










y —m m' J 




{ 


—m m j 



with 



f ' 1 






f 1 1 


Z + l 


\ 














— m 


m 




— m 


rn 


) 





( 1 






( 


1 1 1^ 


+$l>,i-i 






+ h>,i 








—m 


m J 




K 


—m m 
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(A.9) 



Each term in eq. (|A.9[) is given as follows: || 



( I 1 l\ 



\ 



m 



m m 

I 1 Z + l 
— m m j 

I 1 Z -1 

— m m , 



^/Z(Z + l)(2Z + 3)' 



(A.10) 



(-ir 



Z-m+1 



(/ - m + 1)(Z + m + 1) 
^ (/ + l)(2/ + l)(2/ + 3) ' 



l—m 



(-1)'-' 



(Z — m)(Z + m) 
V Z(2Z- 1)(2Z + 1)' 



(A.11) 
(A.12) 



We substitute the 3n — j symbols in eq. ( A.7 ) with eqs. (|A.8|) - (|A.12|) . 
Then we have the form 

F(A,B,C,D: K) 

= J2 J2 J2 J2 K^^fs^u 

{hmi} {hm2} {/3"i3} {Urni} 

x(2/x + 1)(2Z 2 + 1)(2Z 3 + 1)(2Z 4 + 1) 
x {a(Z 4 + 1, m)5 hM+1 + a(Z 4 , m)6 W4 _i} 
x {a(Zi + 1, m)5i 2 j 1+1 + a(Zi, m)5 Wl _i} 
x {a(Z 2 + 1, m)5 /3ii2+ i + a(Z 2 , m)5 fai j a _i} 
x {a(l 3 + 1, m)<5 Z4ii3+ i + a(Z 3 , m), 5z 4 ,z 3 -i}, 



(A.13) 



where 



a(Z,m) = (-1)' 



(2Z - 1)(2Z+ 1)' 



(A.14) 



Because of Kronecker's deltas, the nonvanishing contributions of the 



fourfold summation in eq. (|A.13| ) come from the following six combina- 
tions: 
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q1 


= {h = I, h = I + hh = I + 2, k = I + 1}, 


(A.15) 


ol2 


= {1 + 1,1,1 + 1,1 + 2}, 


(A.16) 


a.3 


= {1 + 2,1 + 1,1,1 + 1}, 


(A.17) 


al 


= {1 + 1,1 + 2,1 + 1,1}, 


(A.18) 


01 


= {1,1 + 1,1,1 + 1}, 


(A.19) 


132 


= {1 + 1,1,1 + 1,1}. 


(A.20) 



Hence, it is sufficient to sum up the terms over the one parameter /. 
The summation of the contributions from the first four combinations, 
I a i,la2,la3,la4:, and that from the last two combinations, lp\, I/32, are 
given as follows, respectively: 



F a {A,B,C,D: K) 



^(2/ + l)(2Z + 3) 2 (2/ + 5) 

l,m 

xa(l + 1, m) 2 a(l + 2, mf 

X \ A l A l+l A l+2 A l+l + A l A l+l A l+2 A l+l 
+ A; A l+1 A l+2 A l+1 + A l A l+1 A l+2 A l+1 



EM 

1=0 



+u l D v l A u l B + ufvfuf 



(A.21) 



and 



F\A, B,C,D:K)=J2 (2/ + lf(2l + 3) 2 a(/ + 1, 



m) 



l,m 
X 



[ A t A F+i A ? a i+i + A t+i A f a< i+i a ? 



U, B U! D + Ul A Ul C 



1=0 



(A.22) 



Here the functions ui and uj, and the coefficient fi and are defined by 
eqs. (p|), ( fnoD , (13) and (H), respectively. 

Summing up ( A.21|) and ( |A.22| ), we arrive at the analytic expression 
of the four-point correlation function in the form eq. (ET 
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Appendix B:. THE THIRD-ORDER NONLINEAR 



SUSCEPTIBILITY OF THE CLASSICAL HEISENBERG 
CHAIN WITH THE PERIODIC BOUNDARY CONDITION 

We write the finite-size scaling function of the third-order nonlinear 
susceptibility for the periodic system as eq. (|2.12|) . We present a deviation 
of this finite-size scaling function in this Appendix. 

We substitute the correlation functions in the right-hand side of eq. 
Q2.2| ) with eqs. (|2.5|) and (|2.7|) . We divide this summation into the sum- 
mation of the four-point correlation functions $ and that of the two-point 
correlation functions \1/ as follows: 



where 



$(T,L) 
¥(T,L) 



1 



3!LT 3 



[$(T,L)-*(T,L)], 



MM 

E(«>r 



(B.l) 

(B.2) 
(B.3) 



Let us show in the following the calculations of $ in details. The 
fourfold summation over the sites ^1,^2,^3 and u in eq. ( |B.2| ) is reduced 
to the threefold summation because of the translational invariance of the 
periodic system. By using the parameters A, B, C and D shown in Fig. 
H], we express the summation of the four-point correlation function, <3>, as 
follows: 

$(T, L) = ±- {MT, L) + 2 (T, L) + 3 (T, L) + 4 (T, L) 



+F(0,0,0,L : K)] 
where Zl is the partition function 



(B.4) 

F is defined by eq. (|A.6|) , and 



h(T, L) = (4-1)! 5l,a+b + c + dF(A,B,C,D: K), 

1<A,B,C,D<L 
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(B.5) 

l 2 (r, L) = 12 2 ^.b+o+u^CO, B,C,D:K), (B.6) 

1<B,C,D<L 



h(T,L) = A E <5 LiC;+Z) F(0,0,C,D:K), 

1<C,D<L 

l 4 (r, L) = sd E ^,c+dF(0, C, 0, D : tf). 

KC.£)<i 



(B.7) 
(B.8) 



Now we calculate each term of ( B.4|) as follows. First, using eq. ( |2.7| ), 
we write 0i explicitly in the form 



L (r,L)= E I 



L,A+B+C+D 



A,B,C,D 
oo 

X 

1=0 



E Af {/, + ^P«P + «? «,V + « W 

1=0 

+ 9l [ufuf + ufuf] } 

E Af {4/M(tf, L) + 2^y 2 (iT, L)} , (B.9) 



where we define V\ and V2 as 

^ g X; 2 Zj ^1 B 1 

^(k,l)^e E E 

A=l 5=1 C=l 

(B.10) 

L-3 L-A-2 L-A-B-l 

v 2 (k,l) = y: e e ^w%w c (B.11) 

A=l B=l C=l 

In obtaining the final expression of ( |B.9| ), we use the invariance of the 
summation of ufvfuf and ufuf under the cyclic permutation among 
A, B, C and D. Each summation in ( |B.1(J| ) and ( B.ll| ) can be carried out 
in the form 



N 

E* r 

n=l 



l-x N )E(x 



(B.12) 



with 
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E{x) 



x 



B 



1 — X 

for | a; | < 1. We thereby obtain V\ as follows: 

L-3 L-A-1 

^i = ^)E E NV 

A=l B=l 



(B.13) 



L-2-A 



|— -E u\-\L - 3) + E^E^) (l - nf" 3 ) 

1 \ 



[£?(«,) - (u l+1 )] £ 



Ul+l 



.L-3 . L-3 



(B.14) 



Here we have used the relation Uj/wj = Ui+i] see the definitions ( |2.6| ) and 
( p.lO| ). In a similar way we have Vi as 



L-3 L-A-2 



^ = %)E E tt^i-uf- 1 - 



A—B 



A=l B=l 

L-3 r 



E 



< (L - 2 - A) + 



1 - 



L-l 



E(^) 2 {(l + ^)(^-3)n 



L-3 



+ (1-u 



L-3 



L-2 



E(ui) 



(B.15) 



The calculations of 02, 03 and 04 are done in the same way. By using 
eq. (|2.7|), we can write down the function 02 in the form 



12 



fc(T, L) = >T Af {ft [V 3 (K, L) + 2V 4 (AT, L) + V 5 (tf, L)] 



+ 9l [V 6 {K,L) + V 5 (K,L))} 



(B.16) 



where we have 



L-2 L-l-C 

\3(A-./.)-E E - u " l '- L - u - r 

C=l B=l 



u t v l u t 



E(u l+ i) [{L-2)uf-ufE(u l+1 ){ui 



L-2 ..L-2 



(B.17) 
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v,(k } l) = j: £ 



C=l B=l 

E( Ul ) (l - ur) - v l E (v?- 2 ~ uf- 



E(v t ) 



L-2 L-l-C 

V 5 (K,L)=Y, E «?«F 

C=l B=l 

= E( Ul ) [E( Ul ) (l - M f- 2 ) — (L — 2)«f" 1 
and 

L-2 L-l-C 
C=l B=l 

— E{ui) [L-2 -£?(«,) (l-^ 2 )" • 
The function 3 is given by 

1 oo 

-</> 3 (T, L) = 2 ]T Af {/, [y 7 (LT, L) + V^LT, L)] + 9l V 8 (K, L)} , 
where 

f 7 (tf, L) = ]T v^~ c = u\E (^f- 1 - tf" 1 ) , 

c=i 

and 

L-l 



C=l 

The function 4 is written as 



1 oo 

-</> 4 (T, L) = £ Af {2/, [(L - + (l - vf" 1 ) 
+ 9l (l + uf) (L-l)}. 
22 



We thereby have the expression of the function $ in ( P-2| ). 



We can calculate the function \1/ in (|B.3|) much simpler. Using eq. 
3J), we obtain the explicit form of \I/ as follows: 



3Zj 
L 2 



.1=0 

oo 



l + U 



L-1 

E 

M=l 



M 



Hi 



U 



.L-M 
I 



£(J + l)Af [l + ttf + 2E( Ul ) (l 



.L-1 



-«=0 



(B.25) 



The above expressions ( |B.2| )-( p.25| ) give a final expression of the third- 
order nonlinear susceptibility. The expressions in the thermodynamic 
limit are given in ref. 3. 

Finally, let us explain the calculation of the scaling limit of the non- 



linear susceptibility, namely eq. (|2.12|) . The scaling limit is the limit 
K — > oo with y = Kj L fixed. It is important to notice that the functions 
ui and vi converge to unity in the scaling limit, except for the terms uf 
and Vi , which turn out to be exponential functions; see eqs. ( |2.13| )-( pTT6| ). 

The functions E{u{) and E(vi), on the other hand, behave as follows 
at low temperatures: 



Ui 



K 



E 



\-ui I + 1 



Ui 

E{vi) 



u, 



Vl 



K 

TTi' 

K 



l-ui 21 + 3 



(B.26) 
(B.27) 
(B.28) 



These terms partly contribute to the prefactor J 3 /T 6 of X3 in Ql -9f) . Since 
we factored out T~ 3 from X3 i n the expression (|B.1|), we pick out the terms 



of the order of J 3 /T 3 = K 3 from the functions $ and (The factor 
1/L in ( |B.1| ) is cancelled out by the factor L in ( |B.4| ).) For example, the 



function V\ in ( |B.14| ) is reduced in the scaling limit to the following form: 



Vi 



K 



K K 



l+l 1 1+2 y 



e -(i+i)/v + 



K K 
I + 1 21 + 3 



++V/y" 
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It is notable that the functions <p2, <p3 and 04 yield terms only of lower 
orders of K than K 3 , and hence vanish in the scaling limit. Thereby we 
arrive at eq. ( f2.12| ). 
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Figure Captions 

1. The spin configuration for the two-point correlation function of 
Heisenberg chain with the periodic boundary condition. Total num- 
ber of spin is L. 

2. The spin configuration for the four-point correlation function of 
Heisenberg chain with the periodic boundary condition. Total num- 
ber of spin is L. 

3. The correction to finite-size scaling of X3T 6 /J 3 for y = J/(TL) = 
1.0, 1.2, 1.4, 1.6, 1.8, 2.0 in the case of the S — 1/2 ferromagnetic 
Heisenberg chain with the periodic boundary condition. We used 
systems with L up to 14. We assumed that the leading correction 
is of the form Jt~/ J. Each solid line indicates the extrapolating 
function based on the four points nearest to the ordinate. 

4. The correction to finite-size scaling of X3T 6 / J 3 for y = J / (TL) = 
1.0, 1.2, 1.4, 1.6, 1.8, 2.0 in the case of the S = 1 ferromagnetic 
Heisenberg chain with the periodic boundary condition. We used 
systems with L up to 10. We assumed that the leading correction 
is of the form JT~/ J. Each solid line indicates the extrapolating 
function based on the four points nearest to the ordinate. 

5. The finite-size scaling function x~3 pen {y) of the ferromagnetic 
Heisenberg chain with the periodic boundary condition, ( |2.12| ). The 
solid curve indicates the scaling function x^ 6 " analytically obtained 
for S = 00. The numerical data, which are extrapolated in Figs. [| 
and £3], are indicated by circles for S = 1/2 and by crosses for S = 1. 
We can obtain the value in the thermodynamic limit, X3 pen (l/ = 0) 
for S = 1/2, by the Bethe-Ansatz method as in (|1 . 1 1|) . 

6. The correction to finite-size scaling of X3T G /J 3 for y = J/(TL) 
= 0.6, 0.8, 1.0, 1.2, 2.0 in the case of the 5 = 1/2 ferromagnetic 
Heisenberg chain with the open boundary condition. We used sys- 
tems with L up to 14. We assumed that the leading correction 
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is of the form JT/J. Each solid line indicates the extrapolating 
function based on the four points nearest to the ordinate. 

7. The correction to finite-size scaling of^T 6 / J 3 for y = J / (TL) = 
0.6, 0.8, 1.0, 1.2, 2.0 in the case of the S = 1 ferromagnetic Heisen- 
berg chain with the open boundary condition. We used systems 
with L up to 10. We assumed that the leading correction is of the 
form Jt~/ J. Each solid line indicates the extrapolating function 
based on the four points nearest to the ordinate. 

8. The finite-size scaling function Xi° pcn {y) °f the ferromagnetic 
Heisenberg chain with the open boundary condition, ( |2.20| ). The 
solid curve indicates the scaling function X3 Pen analytically obtained 
for S = oo. The numerical data, which are extrapolated in Figs. [| 
and|7|, are indicated by circles for S = 1/2 and by crosses for S = 1. 
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M sites 



Fig.1 H.Nakamura, N. Hatano 
& M. Takahashi 




Fig.2. H.Nakamura, N. Hatano 
& M. Takahashi 




Fig.3. H. Nakamura, N. Hatano & M. Takahashi 




Fig.4. H. Nakamura, N. Hatano & M. Takahashi 



Thermodynamic Limit 
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Fig.5. H. Nakamura, N. Hatano & M. Takahashi 



Scaling Limit 




Fig.6. H. Nakamura, N. Hatano & M. Takahashi 



Scaling Limit 




Fig.7. H. Nakamura, N. Hatano & M. Takahashi 



Thermodynamic Limit 
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Fig.8. H. Nakamura, N. Hatano & M. Takahashi 



